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Abstract 

We prove a CLT for skewness and kurtosis of the wavelets coefficients of a sta- 
tionary field on the torus. The results are in the framework of the fixed-domain 
asymptotics, i.e. we refer to observations of a single field which is sampled at 
higher and higher frequencies. We consider also studentized statistics for the case 
of an unknown correlation structure. The results are motivated by the analysis of 
cosmological data or high-frequency financial data sets, with a particular interest 
towards testing for Gaussianity and isotropy. 



1 Introduction 

In recent years, a growing interest has been drawn by infill (or fixed-domain, or high- 
resolution) asymptotics, i.e., the investigation of the limiting behaviour of statistics of a 
single observation of a stochastic processes on a fixed time span or of a random field on a 
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compact space, observed at a greater and greater resolution as the sample size increases 
(see (Stein 1999b) e.g.). 

Such interest was mainly stimulated by various fields of applications. For stochas- 
tic processes defined on subsets of R, the leading motivations have come from finance, 
where it is now customary to observe high-frequency or ultra-high frequency data sets 
collected in a short amount of time (even a single day). In such cases, clearly, the 
standard asymptotic framework, envisaging an ever- increasing time span of observations, 
can be highly misleading and a fixed-domain approach seems a valid alternative. Many 
important contributions have focussed on asymptotic statistical inference on discretely 
observed diffusions, see for instance Kessler and S0rensen (Kessler and S0rensen 1999), 
S0rensen and Uchida (S0rensen and Uchida 2003). 

Our interest here, however, is not on diffusions, but rather on stationary processes, 
as for instance in Stein (Stein 1999a) and (Stein 2002). In this area, strong interest has 
been brought in by the cosmological and astrophysical literature, a particularly active 
field being the analysis of Cosmic Microwave Background radiation. The latter can be 
viewed as a relic of the Big Bang which provides a snapshot of the Universe some 13,7 
billion years ago; as such it is considered a goldmine of information on fundamental 
physics. A widely debated issue relates to the probability law of such radiation: this has 
fuelled considerable activity on testing for Gaussianity on isotropic random fields defined 
on the sphere in three-dimensional Euchdean space. Wavelet analysis has also been 
proposed here, for instance by Vielva et al. (Vielva, Martinez-Gonzalez, Barreiro, Sanz 
and Cayon 2004), Jin et al. (Jin, Starck, Donoho, Aghanim and Forni 2005), Cabella et 
al. (Cabella, Hansen, Marinucci, Pagano and Vittorio 2004). See also (Marinucci 2006b) 
for rigorous results on non Gaussianity testing, but based on the angular bispectrum. 

Wavelet analysis has proved to be a powerful tool in this kind of problems, owing to 
their good localization properties. To the best of our knowledge, however, no theoretical 
analysis exists so far in the literature on the asymptotic (high-resolution) behaviour of 
the random wavelets coefficients for random fields on bounded domains. In this paper, 
we provide some preliminary results in this area. In particular, we focus on isotropic 
random fields on the torus S^; we define a suitable wavelet expansion and we derive 
the correlation structure of the random wavelets coefficients. We use these results to 
establish a central limit theorem for the skewness and kurtosis statistics. These are 
useful to investigate the Gaussianity of the field. As explained before, the asymptotic 
theory is clearly developed in the high-resolution sense. The results are then extended 
to the case where the correlation structure of the field is unknown and estimated from 
the data. 

Our choice of as the domain of the field is a first-step investigation: it is certainly 
desirable to extend our results, for instance, to bounded subsets of M or to S^; such 
extensions are currently under investigation. However also the present results are of 
a practical interest. Indeed, random fields on the circle are the natural environment 
for many geophysical models, for instance concerning atmospheric data (Niu and Tiao 
1995). Also, it is not unusual in the CMB literature to approximate as the union of 
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copies of the so-called ring-torus approach (Wandelt and Hansen 2003). Finally, our 
assumptions can be used for models of stationary time series on a fixed domain. 

The plan of the paper is as follows. In §2 we review the Petrushev construction of 
needlets on general spaces; in §3 and 4 we introduce random fields on the torus and 
the associated random wavelets coefficients; a fundamental bound is established on the 
covariances of the latter. §5 is devoted to a Central Limit Theorem for the Skcwness 
and Kurtosis statistics of these random wavelets coefficients. This result is extended to 
studentized statistics in §6, whereas §7 is devoted to the investigation of the aliasing effect, 
that is, the discretized sampling of the continuous random field of interest. §8 presents 
some Monte Carlo evidence. In the sequel, we use c to denote a positive constant, which 
need not be the same from line to line. 

2 Petrushev and coauthors construction of Needlets 

Frames were introduced in the 1950's (Dufiin and Schaeffer 1952) to represent functions 
via over-complete sets. Frames including tight frames arise naturally in wavelet analysis 
on W''. Tight frames which are very close to orthonormal bases are especially useful in 
signal and image processing. 

We shall see that the following construction has the advantage of being easily com- 
putable and of producing well localized tight frames constructed on a specified orthonor- 
mal basis. 

Definition 2.1 Let M be a Hilbert space, and (e„) a sequence in H; (e„) is a tight frame 
( with constant 1) if : 

V/eH, ||/ir = J]|(/,e.)r 

n 

Let now ^ be a metric space endowed with a finite measure Assume that the following 
decomposition holds 

oo 
1=0 

where the Hi^s are finite dimensional spaces. For the sake of simplicity, we suppose that 
Hq is reduced to the constants. Let Li be the orthogonal projection on Hi : 

V/eL2(^,/x), Li{f){x)^ f f{y)Li{x,y)dfx{y) 

where 

mi 

Li{x,y) = ^e^(x)e^(y) 

i=l 
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mi is the dimension of Hi and (e|)i=i^...^mj an orthonormal basis of Hi. Let us observe 
that we have the following property of the projection operators: 

J Li{x,y)Lmiy,z)dfi{z) = 6i^rnLiix,z) . (2.1) 

The following construction, also inspired by (Prazier, Jawerth and Weiss 1991), is based 
on two fundamental steps: Littlewood-Paley decomposition and discretization, which are 
summarized in the two following subsections. 

2.1 Littlewood -Paley decomposition 

Let be a C°° function supported in |^| < 1, such that 1 > (f){^) > and (f){^) = 1 if 
1^1 < |. Let us define: 

«'(e) = 0(|)-0(O>o 

so that 

V|C|>1, E«'(^) = l- (2-2) 

j 

Actually in the previous sum all middle terms cancel telescopically. Let us define the 
operator 

l>0 

and the associated kernel 

A,(a;, y) = Yl y)= Yl y)- 

l>0 23-l<i<23+l 

Then it holds: 
Proposition 2.2 

J 

V/eH, / = hm Lo(/) + J2 ^•(/) (2-3) 

J— too '—^ 

and, ifMj(x,y) = J2ia(ij)Li{x,y), 

Aj{x,y) = J Mj{x,z)Mj{z,y)dfx{z) . (2.4) 

Proof. 

3=0 3=0 I I 
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□ 



Therefore, as = 1 as soon as 2'^ > I, 

II j]0(^)L,(/)-/ir= Yl \\w)r+ E \\Li{m-<Pi2^ 

I />2J+l 2-'<l<2-'+^ 

1>2J 

(2.4) is a simple consequence of (2.1). 

2.2 Discretization 

Let us define 

m=0 

and let us assume that some additional assumptions are true: 

2. /e J^, gejei^fgejei+i 

3. Quadrature formula: for every / G N there exists a finite subset ^ C ^ and 
positive real numbers > 0, e such that 

V/e^, [ fd^i=J2Kfiv)- (2.6) 



Then the operator Mj defined in Proposition 2.2 is such that Mj(x, z) = Mj{z, x) and 

z I— > Mj{x, z) e J^j+i^i , 

so that 

z 1-^ Mj{x,z)Mj{z,y) e =^3+2-2 , 

and we can write: 

y)^ J Mj{x, z)Mj{z, y) dn{z) = ^ K^^i^. v)Mj{r], y) . 

This imphes: 

Kjf{x) = / A,(x, y)/(y) d^Ji{y) ^ ^ ^)^:;(^, y)f{y) My) 

= E yA;M,(x,77) / y/\,Mj{y,r])f{y)dii{y). 
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This can be summarized in the following way, if we set 

=^2^+2-2 = ^j, i^3,n \/\tMj{x, rj) 

then 

Proposition 2.3 The family is a tight frame 

Proof. As 

/ = Jim (Lo(/) + $]A,(/)) 



2 _ 

J— »oo 



lim ((Lo(/),/) + J](A,(/),/)^ 
■ j<j 

but 

and if -00 is a normalized constant such that {LQ{f),f) — |(/, -^q)!^, then 



□ 



2.3 Localization properties 

Petrushev and coauthors have analysed the previous construction proving that very nice 
localization properties hold. 

In the case of the sphere of R'^"'"^, where the spaces Hi are spanned by spherical har- 
monics, it is proved in (Narcowich, Petrushev and Ward 2006) the following localisation 
property: for any k there exists a constant Ck such that : 

'^^■"^^^l - (l + 2^ arccos(r7,e))^ ' 

In the case of Jacobi polynomials on [—1, 1] with respect to the Jacobi weight, it is 
proved (Petrushev and Xu 2005) the following locahzation property: for any k there 
exist constants C, c such that : 

c2-'/2 

|'0in(cOS^)| < , ■ 

(1 + {2i\e - arccosr/DV^^/slS^cos^) 

where Wo^pin^x) = (1 - x + n-2)"+^/2(l + a; + n-2)^+^/2, -1 < x < 1 if a > -i, P > -\. 
In the following section, we consider the case, which is our framework, where ^ is the 
torus, (ejfc) is the Fourier basis, and if^ = Span{e^}. 
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2.4 Quadrature Formula on the torus 
Proposition 2.4 Assume that ^ — T, the torus. If, for m eN, 



\k\<m 

the quadrature formula (2.6), holds for 



1^^, I e {0,...,m}|, 
I m + 1 J m + 1 



Proof. Let T be the torus, identified with [0, 27r] and endowed with the measure f dfi = 
'h lo^ •^(•^) / : T — i> C with the following expansion on the trigonometric basis, 

-f \ ^ „ Akx 

f ^ J-^dke . 

k 

For m e N let us define, 

^ m 

It is obvious that T^+if is periodic with period Therefore T^+i/ has the expansion 

— / T^+if{x) e-'^'' = if A; ^ O(mod m + 1) , 
i- ^ T^+i (/) (x) e-^'^^^+i)- = afe(^+i) . 

Hence, 

If / is a trigonometric polynomial of degree smaller than or equal to m, we have 

m „ 

Therefore, if / G J(^m^ 



1 1 /■ 



m , _ 

i=0 



□ 
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2.5 Localization properties for trigonometric series 

Following the steps of the sections above, we have for 77 e =^j+2_i, that is 77 = k e 
{0,..., 2^+^-1}, 

= (2-7) 

We prove now at the same time the concentration of the wavelet and of A^. Let be a 



-2 



— TT 







TT 



Figure 1: Comparison between the wavelet ipjQ (solid) and the trigonometric function at 
comparable frequency (dots). Here j — 4. The local character of the wavelet is apparent. 
In this picture the wavelet has been normalized in order to have the same norm 
(otherwise it would be smaller). 

a continuous compactly supported function and let us consider, associated to A and j, 
the function 

e,(x) = 5^A(X)e^'^. 

We denote by — W^{W) the Sobolev space of functions with integrable weak deriva- 
tive of order k. 

Theorem 2.5 Let A be a continuous, compactly supported function such that A e 
for some k >2. Then, for a// m e N there exists a constant Cm,k such that 



r 1,2 ('"+^)-' 
I sjv J\ - (1 + |2J2;|)'= 



(2.8) 
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Proof. Clearly ^j{x) is a trigonometric polynomial. Moreover let us put 

B{x) = ^{A){x) = i- 1 A{y)e'y^dy. 
If e Li(M), by the Poisson summation formula : 

k k lez 

and more generally if e Li(R) 

D'^iQix) - 27r^2(™+i)^D™(5)(2^(x - 27r/)) . 

But 

D"^{B)ix) = D^.^{A){x) = ^i{tyrA{y)){x) 

is bounded (and, by the Lebesgue-Riemann lemma, even vanishes at infinity). Further- 
more, 

{-ixfD"^{B){x) = {-ix)'^^{{iyrA{y)){x) = = 

k 

and this function is bounded as ^4 e W^. Therefore 

( 



Hence ^ 

The result is now a consequence of the following lemma. 
Lemma 2.6 Fork>2 

is a 27r-periodic function such that 

ej{x) < ^ 



(1 + |2^x|)* 
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Proof. Let |i| < TT. Then 



5(1 



{l + \2^{x-27rl)\)'' (l + \2jx\) 



1 1 

1 + \2^x\)'' ^ (1 + \2^{x - 27r0|)^ 



Since 



Sit 



1^0 



+ \2^{x-27rl)\y 



< 



+ 2^(2/- l)7r; 



< 



;i + 2%)fc 



1 + 2% 



< 



(1 + 2i7r) 



1 + 



Z>2 

°° M + 2% 



1 + 2^X71 



+ 2i(2/- l)7r 

fc 

dx 1 < 



< 



(1 + 2%)'= 



one gets finally 



14 5 



□ 



Remcirk 2.7 We are well aware that other choices of a well localized basis is possible. 
As for instance the Lemarie- Meyer wavelet system (Meyer 1992), that is also formed by 
trigonometric polynomials and that is moreover ortho normal, whereas frames are not 
orthonormal and redundant. Recall however that, as it is well known, the advantage of 
redundant frames is robustness. Actually, assume that {ei)i^i is an orthonormal basis of 
an Hilbert space H, and that we have a noisy observation of a function / e by 



(2.9) 



where (ej)i is a sequence of r.v.'s modeling the noise. Then the error of the estimator 
X^igj y^ej is easily computed 



iei iei iei 

Conversely, if (ei)jg/ was a tight frame then it is easy to check the inequality 

2 



^|(0,e,)r = ||0ir= <J2\\ 



iei 



Therefore the error is smaller 

iei iei 
In fact often the inequality appears to be strict. 
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3 Assumptions and random wavelet coefficients 

3.1 Assumptions on the model 

Consider the random field 

oo 
(=— oo 

where 

oo 

= 0, Ewi = 0, E\wi\^ ^Cu 1^1,2,...,^ Ci< oo . 

l = —OD 

Throughout this paper, we introduce the following regularity condition on the behaviour 
of the angular power spectrum. 

Assumption Al There exists a function g : N ^ [ci,C2] such that g e for some 
M > and 

d = gil)!"^ for alU e N , a > 1 . 
For some results to follow, this assumption is strengthened to 

Assumption A2 Al holds and there exists a sequence of functions hN{u) : [|,2] — > 
such that 



1-C2 ' ClJ 



M#):=^, ^</<^,7V = 8,16,32, 



2 

^ \ 4 / 

and 



sup sup < C , some C > 0, some M e N , (3.10) 

^ 1/2<M<2 

where h^^^ denotes the M-th order weak derivative oi h^- 

Remcirk 3.1 Condition Al is a very mild regularity condition; it is imphed, for instance, 
if g{l) is any trigonometric polynomial bounded away from zero. The requirement a > 1 
is necessary to ensure the sequence Ci to be summable, which in turn is a consequence of 
the finite variance of the field. Condition A2 is a slightly stronger smoothness condition, 
which implies that /ijv G Wf^. 

3.2 Random wavelet coefficients 

We recall the frame introduced in (2.7), namely 

^..(^) = ^ E «(^)^'^^-''^ v = ^,,ke{Q,..., 2^>^ - 1}. 

1^0 
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The notations will be shortened into the following way. For j e N, we put N — 2-^"'"^, 

27r 

i^Nkit) = ^N{t -kr), T=—,ke{0,...,N-l} 



oo 



t=— oo 



where a is a "^"^ function, compactly supported in [|,2]. Hence, ^n{1) has support in 
(f,f); indeed 

^N{l) = ^a{^). 

We have also 

M) - ^ f_ V'Jv(^ + kT)e''^ = e-'^'^i^Mil) . 
We now define the associated wavelets coefficients of the process X 

^Nk-^^j X{d)i^^k{'d)d'd^^ j X{'d)i^N{'d + kT)d'd, A; = 0,1,..., TV -1. 

It is immediate to see that EI3Nk — 0; also 

E{l3Nk,l3Nk,) = f r EX{'d)X{'d')^^k,WNkM)d^d^' - 

—IT J —IT 

EX{^)X{^')xljNk,mNk,{^') d-d d-d' = 

OO OO 

= E Cii,^kMNk,{-l) = Yl Ci\ip^{l)\'e'^^^'^-'^^ = 



(3.11) 



1 aa(#)^e^M.-'=^). 



AT 

Z=— oo 



Next we study the asymptotics of the correlation coefficient, defined as 

^N/8<\l\<N/2 ^l^'^il^) 



27r 

n m Q _ ^N/8<\l\<N/2'-'l(^ Kn)^ 

Uorr [PNki:PNk2) — ^ 2(41^ 



^N/8<\l\<N/2 ^I^'^^n) 

Lemma 3.2 Under Assumption A2, 

|Corr(/3,.,„/3^,,J| < ^ ^kT-k,]^/,^ 
for some cm > where [a\b means a(mod)&. 



3. Assumptions and random wavelet coefRcients 



13 



Proof. Write 



Corr(/3jvfei,/3jvfe2) = 



27r 



^ I]Ar/8<|/|<Ar/2 Cm/ a"' ^ N 



N ^N/8<\l\<N/2 Cn/4"' Vat 

and note that under Assumption A2 it holds for the denominator 

C, 



< ^ E 7^«'(#) ^ ^22" , as ^ oo 



N/8<\l\<N/2 



Thus we focus on 



- E 



AT ^ Cjv/4 



27r 



< 



< 



,27r, 



^(0 (Al\-a„2(il\ ilj^{ki-k2) 



Ar/8<|i|<iV/2 

To complete the argument, we use Theorem 2.5. To apply the theorem, we notice that 
for all N 

An{0 := /iiv(0(0"V(0 

is the (sampling of) the Fourier transform of an infinitely differentiable expression, so we 
obtain 



N/8<\l\<N/2 

1 2ir 

\^ h ('4^W4^^-a_2/4^^^^^-i^7(fel-fe2) 



,2-K , 



N 



< 



N/S<\l\<N/2 

2cmN 



< 



< 



2c 



M 



AT (1 + Ar(f [k - k%^,)r " (1 + - k'U/,r ' 



which gives the required bound; note that cm does not depend on A^, in view of (3.10). 

□ 

Remcirk 3.3 Lemma 3.2 highhghts a quite remarkable property of random wavelet co- 
efficients. Indeed it entails that wavelet coefficients located at finite distance are asymp- 
totically (with respect to the frequency A" = 2-^+^) uncorrelated. 



We write 



N 

N/8<\l\<N/2 
^Nk 



N/S<l<N/2 



(3.12) 
(3.13) 



so that -B/S^fe 
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Remcirk 3.4 It holds, as A?" — > oo, 



N/8<l<N/2 



dt 

1/2 



therefore, as under Assumption A2 < ci < Ci/CN/i < C2 < +oo, there exist constants 
< < such that 



•-1 ^ /-( — '-2 



for every A^ > 0. 



In view of the asymptotic results of next section, we shall always focus on the Gaussian 
motivated by our statistical applications. 

Assumption B The field X is Gaussian. 



4 Asymptotics of the wavelet statistics 
4.1 The sample mean 

Our first aim in this Section is to investigate the asymptotic behaviour of the sample 
mean for the random wavelet coefficients. More precisely, define 

1 ^ ^ 

Mn -j^^f^Nk ; 
fe=i 

we have immediately EMn = 0. Under Assumption Al, it is also simple to show that 

1 ^ ^ ^ 



N/S<\l\<N/2 " k=l 



'^''^=' ^ (4.14) 



•Ml. 



2 

= 0, 



from the well-known properties of the Dirichlet kernel 

.27r/, 

2J iv = for aU / e N such that ^ ^ 2k7r , k = 0, ±1, ±2, . . . . 



k=i 



It follows immediately that Mjv = with probability one. It is interesting to realize 
what happens here. Given the fast decay of the covariances established in the previous 
Section, we might have expected standard asymptotics to go through, i.e. a Central Limit 
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Theorem for the normahzed sample mean. This is not the case because the variance is 
degenerate; intuitively, this is due to the wavelet transform which is 'overdifferencing' the 
random field. Put in another way, if we view the wavelet coefficients as a discrete time 
periodic random process, then this process as a zero spectral density at the origin. This 
complicated dependence structure does not prevent, however, the Central Limit Theorem 
to hold for higher-order statistics, as we shall show in the sequel. 



4.2 Skewness and Kurtosis 

Motivated by testing for non-Gaussianity on random fields on the torus, we introduce 
here the Skewness and Kurtosis statistics of the wavelet coefficients. More precisely, we 
consider (recall (3.12) and (3.13)) 

1 ^ ^ 1 ^ ^ 

^^■■--7^T.^Nk and [/^ ^(/?^, - 3) . 

^ k=l ^ k=l 

We have immediately 

E[Sn] = E[Un] = . 

The variance of these statistics is more complicated; in view of the following Lemma, it 
is convenient to introduce Fejer's kernel, defined by 



knU) := 



1 sin^(iA^f) 



27rN\^ I 27rN sin' at) 

k=l ' 

Fejer's kernel vanishes at the Fourier frequencies ^ /, unless / = kN , A; e Z; moreover, if 
/ = kN, Kn{^1) = §. We define also 

4. := ^4 E Q,a\f)Q,a\f)Q,a\f)K^{{k + k + k)r), (4.15) 
^ hhh 



and al^ = aj^^ + a^^^, where 

-■=^W^I1 Q,a\f)Q,a\f)Kr,{{h + W , (4.16) 

^ hl2 

'^^ hhhu (4-17) 

XKN{{ll + l2 + h + k)T). 

Remeirk 4.1 As Fejer's kernel vanishes at the Fourier frequencies ^l, I kN , A; e Z, in 
the previous expressions most terms vanish. Taking into account the fact that a^{x) — 
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unless ^ < \x\ < 2, one can write in a more computationally tractable way 

hh 



Likewise 



2 '''2 1 V^^2^4/4Zi\ 

'^iUn ■- 'ZA-f^Z^'-'h^^ KIT)- 



A similar, a bit more complicated, expression can easily be derived also for cr2UN- 
Lemma 4.2 Under Assumptions Al and B 

E[Sl] = 4. , (4.18) 
E[Ul] = <^lu.+<ylu. ■ (4.19) 

Proof. For (4.18), we note that, by the diagram formula (see (Adler 1981), p. 108 e.g.) 

1 ^ 

^ fcl,fc2 = l 

1 ^ 

= -T^— 6" E {S-^'f/^iVifci/^ArfeJ-fi^f/^Arfea/^Arfeal-^'f/^Ariki/^Arika] + 

Q iV R ^ 1 3 

-'^ fci,A:2=l " fci,A:2=l Z 

We have seen already in (4.14) that the first term is equal to zero. As for the second one, 
we obtain (recall that t = ^) 

N 



^ fcl,fc2 = l 



N 

AWT E {4 E a.a^(f)a<.^(f)C„a^(f 

^ fci,fe2=l hhh 

127r 1 



^ hhh 
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whence (4.18) follows. Likewise, for (4.19), we have 

E[0Nk, - 3)(&. - 3)] = E0j,J%,J - 9. 
Again by the diagram formula and recaUing that E[/3jf,^_^] = 1, 



so that 



Now 



1 ^ ^ ^ 

fel,fc2=l 

24 ^ 4 72 ^ 

^ fci,fc2=l ^ ki,k2=^ 



1 1 4 



^ fcl,A:2=l ^ fci,fe2=l 



1 ^ 



fei,fc2=l lihhU 

E Qi«'(f )a2a'(f )Q3a'(f )Q4«'(f )i^iv((/i + /2 + /a + /4)r) 



7V4 



On the other hand 

AT N 



^ fcl,fe2=l ^ fcl,fc2=l « 



27r 



W2 

□ 

Remcirk 4.3 By the same arguments of Remark 3.4 it is immediate that, as A" — > oo, 

^ E )^iv(t (^1 + ^2 + /3)) - 



N3 

a^(ti)a^(t2)a^(t3) . ^ ^ ^ dtidtidta . 
1/2 7i/2 Ji/2 2 sm [ti +t2 + ts) 



Therefore under Assumptions Al and B we have 



< ci < (7|^ < C2 (4.20) 

for some constants < Ci < c^- By the same argument we see that also afj^ is bounded 
and bounded away from 0, for every N > 
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5 The Central Limit Theorem 

This Section is devoted to the Central Limit Theorem for our statistics of interest. The 
idea is to prove the results by the method of moments. To analyze the behaviour of 
higher order moments, we shall extensively use the already mentioned diagram formula, 
that states that, for a multivariate zero-mean Gaussian vector {Xi, . . . , X2k), it holds 

£;(XiX2...X2,) = J]E(X,,X,J...E(X,,,_,X,,J . (5.21) 

Consider the cartesian product / x J, where J, J are sets of positive integers of cardinality 
7^(/) = P, i^{J) = Q; it is convenient to visualize these elements in a PxQ matrix with P 
rows. A diagram 7 is any partition of the P ■ Q elements into pairs like {(ii,ii), (^2,^2)}; 
these pairs are called the edges of the diagram. We label r(/, J) the family of these 
diagrams. It can be checked that, for given I, J, there exist (P-Q — 1)!! different diagrams, 
each of them composed by ^P-Q pairs; we recall that {2p — 1)!! := (2p — 1) • (2p — 3) • • • 1 
for p = 1,2,... 

To any diagram we can associate a graph having / as the set of vertices (or nodes) 
and connecting any two of these vertices, ik,ik', by an arc whenever in the diagram an 
edge of the type {{ik,jk), {ik'jjk')} is present. This graph is not directed, that is, (^1,^2) 
and {i2,h) identify the same arc; however, we do allow for repetitions of edges (two rows 
may be linked twice). We shall use some result on graphs below; with a slight abuse of 
notation, we denote the graph 7 with the same letter as the corresponding diagram. 

We say that 

a) A diagram has a flat edge if there is at least a pair {(ii, ji), (^2,^2)} with ii — 12. 
we write 7 G Tp{I,J) for a diagram with at least a flat edge, and 7 G r;p(/, J) otherwise. 
A graph corresponding to a diagram with a flat edge includes an edge ikik which arrives 
in the same vertex where it started. 

b) A diagram 7 G r-p(/, J) is connected if it is not possible to partition the corre- 
sponding graph into two sets A, B such that there are no arcs connecting the nodes in A 
with nodes in B. We write 7 G Tc{I, J) for connected diagrams, 7 G r^(/, J) otherwise. 

c) A diagram 7 G J) is paired if, given any two sets of edges {(^i, Ji), (^2,^2)} 
and {(^3,^3), (^4,^4)}, then ii = i^ implies 12 = i^; in words, the rows are completely 
coupled two by two. We write 7 G Tp{I, J) for paired diagrams. 

Obviously if P > 2 a paired diagram cannot be connected. Note that if Q is odd, 
paired diagrams cannot have flat edges, so that the assumption 7 G ^-p{I, J) becomes 
redundant. If / or J (or both) can be simply taken as the set of the first p or q natural 
numbers, i.e. / = {1, . . . ,p}, J = {1, . . . , g} we shall occasionally write r(/, q), Tp, J) or 
r(p, q) for r(/, J). For a nice and comprehensive discussion on the diagram formula, we 
refer to (Surgailis 2003); see also (Nualart and Peccati 2005) for recent advances in this 
area. 
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Theorem 5.1 Under Assumptions Al and B, as N ^ oo, 



Hence 




ifpi,P2 are even 
otherwise. 



Proof. For brevity's sake and notational simplicity, we focus on the case p2 — 0; the 
general argument can be pursued under the same lines. The idea is to use (5.21) and 
partition the various summands in this expression according to the nature of the asso- 
ciated diagrams/graphs. To this aim, wc visualize our coefficients as positioned on a 
diagram with / = 2pi = 2p rows and J = 3 columns; we obtain 

N ki,...,k2p 

E M n {^^ E 

'^N ki,...,k2p ke{ki,...,k2p} N/8<\l\<N/2 

= e[ n n E ».«(fy"] = 

"iV kx,...Mp fce{fci,...,fc2p}i=l Af/8<|«|<Af/2 

= 1^ E E n { E «(#)<.(f)e*'"e"'--} 

fci,..,fc2p7er(2p,3){(mj)(m',j')}e7 N/%<\l\,\l\'<N/2 

= E E n { E 

'^N fci,..,fe2p7er(2p,3){(m,i)(m',j')}e7 N/8<\l\<N/2 

where wc used the fact that ElwiWii] = unless I' — —I. For fixed ki, . . . ,k2p, let us 
concentrate on the contribution of a single diagram 7 e r(2p, 3). 

n { I] Cia\§)e''^''^e-''^''^''j = 

{{m,j)(m',j')}e7 JV/8<|/|<Ar/2 
'll,---,il3, — ,'2p,lr--,'2p,3 {{m,j){m' 

where the indices Imj, m = 1, . . . , 2p, j = 1, 2, 3 vary between and N/2. Now, as 
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every km appears exactly three times in the diagram, 

{(m,j)(m'j')}e7 



{(mj)(m'j')}e7 {(fnj){"j'J')}e7 

2p 



{(mj)(m'j')}67 "*=1 

where l-mj-a = hn,j if {nT-:j) is a departing point in 7, /mjj-y = — ^mj if is an arrival 

point. Summing on the possible values of /ci, . . . , k2p, we get 

/ J J. J. ^ N > 'mjj'm'j' 

fei,...,fc2p {(m,j)(m'J')}e7 

n <^'^.>'.'./) ( n ^iv([/m,l;7 + ^-,2;^ + im,3;7]T)) = 

{(mj)(m' ,jr')}e7 "»=1 

X 



2p 



{(m,i)(m',j')}G7 {(mJ)(n*' j')}e7 

2p 



='5(7;'i,i>---i'2p,3) 

m=l 

= 5(7; . . . , /2p,3) n (^^^.«'(%^)^^™.«'(%^)^^^,3«'(%^)) X 
m=l 

Xl?jv([^m,l;7 + ^m,2;7 + ^m,3;7]'?'). 

Let us define 

= jy2^3 (^'m,i«^(%^)C'im,2^^(%^)^'m,3«^(%^)) -^^AT ( [^m,l;7 + W;7 + W;7]'^) ■ 

In conclusion 

2p 

^^N = ^ ^(7; ^1,1, • • • , ^2p,3) n ^'m,l,«m,2,«m,3;7- 

7er(2p,3) /i,i,...,«2p,3 m=l 

Recall that 5(7; /i^i, . . . , /2p,3) = unless lm,j = Im'j' for every {m, j){m' , j')} G 7. The 
proof is done by proving that the leading contribution to ES'^ is given by paired diagrams 
whereas the non paired ones are negligible in the asymptotics. This is made explicit in 
the two following lemmas. 
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Lemma 5.2 For the terms corresponding to the paired diagrams 7 e rp(2p, 3) we have 



Proof. Remark first that the number of possible ways of partitioning the 2p rows of 
the diagram into subsets of cardinahty 2 is exactly (2p — 1)!!. Also, in every diagram 
7 e rp(2p, 3) the contribution of two paired rows is exactly ag^. 



To conclude the proof, we need only to show that the terms corresponding to all remaining 
diagrams 7 ^ rp(2p, 3) are asymptotically of smaller order. 

Lemma 5.3 All terms corresponding to diagrams with connected components of order 
larger than ^ (7 ^ Tp{2p, 3)) are of order 0{^^^). 

Proof. We focus on any two nodes that are connected but not paired; such two nodes 
certainly exist, otherwise 7 G Fp. We consider the case where there is a single edge 
linking these two nodes; the proof in the remaining case is entirely analogous. Without 
loss of generahty, we label edges and vertices in such a way that the edge connecting 
these two nodes is labelled {(1, 1), (2, 1)}. 

As in Lemma 3.1 of (Marinucci 2006a), we apply iteratively the Cauchy-Schwartz 
inequality to show that 



2p 



5] ^ Sir,h,u...,i2p,s)ll^i'. 





□ 




2p 



1/2 



where 




'1,1,(2,1 



On 

X — ^ Cil,ia^(%^)-Div((^l,l;7 + ^l,2;7 + ^1,3;7)t)-Cjv((-/i,1;7 + ^2,2;7 + ^2,3;7)'r) ■ 
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Now 



/ V lm,llm,2lm,3;7 
llm,2l m,3 



m,3 



o 



,1;7 ^" ^m,2;7 ^" ^m,3;7 

]r) = 0(1) . 



On the other hand 



y 



V y2 

/ V il,2il,3'2,2(2,3;7 

h, 2^1, 3^2, 2^2, 3 

ChoO^i %^ ) C;, 3 ( %^ )Cu^a'^{ %^ ) C/2 3 a 



X 



'N 



'l, 2(1,3^2, 2^2, 3 

J_ >sp DN{{h,l;^ + ^1,2;7 + /l,3;7)^)-DjV ((-^1,1;7 + ^2,2;7 + ^2,3;7)t) 



X 



'1,1 



'AT 



N 



Now we observe that 

-DAr((/l,l;7 + + /l,3;7)T)DAr(( — + /2,2;7 + ^2,3;7)t) 



N N 

— ^ ^ g('l,li7+h,2i7+il,3;7)Tfel ^ ^ g(— /l,l;7+i2,2;7+i2,3;7)''"fe2 

A;i=l k2=l 
N N 

— ^ ^ g(il,l;7+il,2;7+h,3;7)'rfel ^ ^ g(il,l;7+i2,2;7+i2,3;7)'r(fel -tl) 
fel = l U=l 

Dn{{Ii,2;j + ll,3;j + ^2,2;7 + ^2,3;7)''') -Dat (( — ^1,1;7 + i2,2;7 + ^2,3;7)'''), 
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whence 

«^ ^ ^ Al 



< 



r-DAr((/i,2;7 + /l,3;7 + ^2,2;7 + ^2,3;7)''" ) X 



1 

TV 

^ hv ^ ZT^hl;^0''^{^)DN{{-h,l;^ + ^2,2;7 + ^2,3;7)'r ) 



< 2 Iv ^ < 

Kl,2;7 + ^1,3;7 + ^2,2;7 + ^2,3;7| + I N f-^ I — + /2,2;7 + ^2,3;7| + 1 



< 



Kl,2;7 + ^1,3;7 + ^2,2;7 + ^2,3;7l + - , , I <'i,i;7 ' •':^,2;7 

''1,1;7 

<c, 

Kl,2;7 + ^1,3;7 + ^2,2;7 + ^2,3;7| + 1 



Thus we can conclude that 

1/2 

,3^2,2i2,3;7 



( ^^i,2'i,: 

h,2il,3i2,2i2,3 



< 



^ ^ 2h'k.^.,3 (Kl,2;7 + ^1,3;7 + /2,2;7 + /2,3;7l + 1)7 



□ 



6 Studentized statistics 
6.1 Estimation of g\. 

The statistics described in the previous sections can be impossible to compute in practice, 
as they depend on the correlation structure of the field, which is in general unknown. In 
this Section, we show how asymptotic variances can be consistently estimated from the 
data, in the presence of observations at higher and higher resolution. We start from the 
variance of the wavelets coefficients, which we recall is given by 

N/?,<l<N/2 

An obvious estimator is provided by 

N/S<l<N/2 

Of course, a% is unbiased and mean square consistent for cr^, in the trivial sense that 
both converge to zero as N diverges. The following result is stronger. 
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Lemma 6.1 Under assumptions Al and B, as N ^ oo, we have 

lim ^ = 1 m L2 . 
Proof. It is immediate to see that 



and 



N/8<l<N/2 ^t"^ 

Var(a^)= J] Var(|^,| V(#) = 2 J] Cfa\f). 

N/8<l<N/2 N/8<l<N/2 

It is sufficient now to note that, under Assumption Al, 

< ci < < C2 < oo, for N = 2^+2 and for all I e [A^/8, A^/2]) , 

Cn/4 

whence 



Var 



N' 



Remark now that ~ ^ jl^^a^{t)dt, E/«^(f) ~ ^ /i/2«^W^^> so that the 

□ 



left-hand term tends to as 



6.2 Estimation of the variance for Skewness and kurtosis 

We now go on with the estimation for the sample variance for the statistics Sn and 
Un- We note first that under Gaussianity {\wi'f /Ci)i is a sequence of independent and 
identically distributed exponential random variables with mean unity; we define for all 
p e N 



i=h 



\ 7—;. 

to be quite explicit we have for instance 



\Wi\ 

C, 



|2 





if 


l^ll ^ 1^2! 




if 






if 




I 


if 


among \li 




if 


— 1^2! 
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In view of (4.15)-(4.17), a natural proposal is to consider for the Skewness 



127r ^ 1 



and for the Kurtosis 

^k-^iu^+^2U^ ' (6-23) 

^1^. ^= 5§E;^l^'.l'«'(f + ' (6-24) 

(Tat ^ i^i^ 

^"^^^■■=^% E ^|nkdV(#)}x..(t(Z, + /. + Z3 + Z4)). (6.25) 

hhhh 'i'2/3/4 L i^i^ J 



wi,\'aYj^)\wifa\f)\wifa\f)Kr,C^{h + I2 + k)) (6.22) 



'N 



RemcLrk 6.2 Using the properties of Fejer's kernel recalled in §4.2, in the summations 
above most terms vanish. From a computational point of view more tractable expressions 
can be derived in the spirit of §4.2. In particular it holds 

Lemma 6.3 Under Assumptions Al and B, as N ^ 00, we have 

P ^1U% P ^2U% p 

^ ' 2 ^ ' 2 ^ ' 

Proof. We give the proof for S'^j;^ , S'lj;^ only, as the remaining case is entirely analogous 
(indeed slightly simpler). Let us denote ^fu^^ — o'iun ' ^n/^nj ^^^^ same as in 

(6.22) with CTjv replaced by a^. As E (^j^\wij^\^\wi^\'^^ — Ci^Ci^ for every Zi, ^2, it is clear 
that 

sf^^^l. (6.27) 



Moreover, by the alternate expression (6.26) and in view of Remark 3.4, 

As we know that under Assumption Al (yluj^ is bounded away from zero (see Remark 
4.3), this implies that Var(a^^^/(T^^^) ^ as — 00. 

The argument for o^^jj^ is similar; indeed if we define ^2Un analogy with ^lu^^ , it 
is immediate that E[a2Uj^/cr2Uj^] = 1. On the other hand, note that the summands in 
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(72^^ have a martingale-difference structure on the lattice (see Poghosyan and Roelly 
(Poghosyan and Roelly 1998) e.g.), whence, in view of Assumption A 



Var 



/ ^2 



1=1 



KNifih + k + h + k)) 



N 

k 



-^iw^H ^^(tai + ^2 + /3 + ^4))Var|n^| 
^ hhhU ^ l=li ^ ^ 

=o(]^Ev4n^})=o(^)=o(i). 

IMsh ^l=h 

We have thus proved that o^fuff/'^iUN ~^ ^ ^IunI^^Un ~^ 1 as — > cxd in L^. 
Therefore, in view of Lemma 6.1, o^Iu^/^^iun ~^ ^ ^2Un/^2Un 1 a-s A" — > oo in 
probability. The rest of the proof is quite similar. 

□ 

As an immediate consequence of Theorem 5.1 and Lemma 6.3 we have the following. 
Theorem 6.4 Under the assumptions Al and B, as A" —> oo 



-' N I N-^oo 



N{0,l2) 



7 Aliasing 

The tests provided in the sections above are based on the wavelet coefficients 

PNk - ^ / X{^)iljNk{^)d^, . A; = 0, 1, . . . , AT-l 
In practice, P^k will be approximated by the sums 

^ M-l 



M 

m=Q 

The purpose of this section is to prove that, if our data have enough high frequencies, 
this approximation does not affect the test. 

We need to strengthen our previous assumptions as follows. 

Assumption C As Z ^ oo we have 

Ci = L(or°, a > 1 
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where L{1) denotes a slowly varying function (Bingham, Goldie and Teugels 1989), which 
we assume to be bounded and bounded away from zero. 

Assumption D A'" is such that 

T7 + /, ^ as M ^ OO . 



Assumption C is mild, entailing simply a regular behaviour of the angular power spectrum 
at infinity. Assumption D is a sort of bandwidth condition, suggesting that the frequen- 
cies that we can use fruitfully for statistical inference must grow more slowly than the 
samphng rate. The condition become less and less tight the faster the angular correlation 
function decays to zero: for instance if a = 4 we must impose N = o(M^/^). In practice 
a can be estimated from the data; the most cautious choice can be M = o{-\/N), as q; > 1 
is implied by the finite variance of the field. 
We have the following result. 



Proposition 7.1 Under assumptions^Al,B,C and D the result of Theorem 5.1 remains 
true when replacing the 13 nu 's by the Pnu 's. 

Proof. It holds. 



K|>M-f 



m=0 



Therefore, 

E\f3Nk — PNk\ 



E 



^ M-1 



'\>M-i 



m=0 



2-Km 



where 



1 

M2 



M-1 



m=0 



< 



+ W 



— for all /c > , some > , 



in view of Theorem (2.5), which implies 



M-1 



m=0 



< 



\l + M[l/M]2^^ 
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Under Assumptions C and D we have easily 



Ci\Am(i)\^ = o(cm ^2 7r~l^^(0r] 

\e\>M-f ^ \e\>M-f ^ ^ 



oo M ^ 



u=\ v=—M 

oo M 



M 



M{uM+v) I 



^ u=l v=-M ^ II// 



and 



i^(|/37Vfc-/3. 



Nk\ 



E\l3Mk\ 



<c^ <c( — ] 



(7.28) 



The result then follows from the following lemma. 

Lemma 7.2 For and Y^^ mutually Gaussian centered random variables, set 



Let us assume that the sequence 



'Var(Er=i^l) 
converges in distribution to a variable X, where for 71, 72 > 



(7.29) 



71 



< 



If"! 1 

^ -^^^1 E^i*"! < 72 and Cn = o(-) . 



Then 



(7.30) 



(7.31) 



a/so converges in distribution to X . The same result is true if we replace in (7. 29), (7. 30) 
XI by XI - EXl and replace m (7.31), Yl by Y^^ - EY,l 

Proof. We shall actually prove a stronger result, namely 

/-^i=l in Z-^i=l in 



lim E 



= . 
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Using {x — yY — — — 3x{x — yY + 'ix^{x — y), we get : 



1 1 ,4,1 



< E\X,^ - YJ' + 3[E\X,^ - Y,X]^E\X,r,\']^ + 3[E(X,, - Y,^y]^E\X, 
Now, when Z is a Gaussian random variable, for h e N*, 



2 



where o"^ is the /i-moment of the standard gaussian distribution, centered and with vari- 

h ^ 

ance 1. Therefore, we have E\Xin — Yin\'^ < ahCn [EXf^]^ , and 



mi 

n 



2]^ 
ini 



i=l i=l i=l i=l 

n 

< J2[EXUH<^3cl + 3a|c| + 3a|c„}. 

1=1 

Therefore 

EkJ^^M < ^Skt^^lM < = 0(1), as n ^ 00. 

This proves the result for the third power. As for the forth one, we proceed in the same 
way, and prove using the same path, 

— y'^l < 4|x|^|x — y\ + 6x'^\x — y\'^ + 4|a;]|a; — y\^ + \x — y\'^ 
E\Xt-Y,i\<C^[EXl]' 

En rp I Y'^ — I 
i=l'^\^in ^ir 



< Cy/C^^ 0(1) . 



□ 

Our final result in this Section extends the analysis of the aliasing effect to studcntizcd 
statistics. In particular, in the previous section the variances of Skewness and Kurtosis 
where estimated on the basis of the spectral coefficients {wi} , which are obtained as 

wi^^ r X{'&)e-''U^ , Z = 1,2,...,7V. 

As before, in practice these Fourier coefficients will be approximated by interpolations 
sums, i.e. for M > N we have to consider 



M-l 



. 27rm , 



m=0 
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We note that 

M-l oo , M-1 oo 



m=0 fc=l m=0 k=l 

^ oo oo 



(7.32) 



M 

fc=i fc=i 



It follows, using Assumption D, that 

C, = 2. ^^+fcM < C-— /c <c(^-j =o(^). (7.33) 

' ' k=l * ifc=l 

Proposition 7.3 Under Assumptions Al,B, C and D the result of Theorem 6.4 remains 
true when replacing the wi 's by the wi 's. 

Proof. It is clearly enough to show that 

hm - hm - (7.34) 

in probability, where ^"sj^i^Un defined in (6.22) -(6.23), 

and = alu^ + where 



' - ^^"^ E7^l^'il'«'(f)l^'^l'«'(t)l^'3|V(f)X;v(t(il + /2 + i3)) 
hhh ^1^^'^ 



-^|^E7^i-'^i^«'(t)i-^.|V(t)i^iv(tai+^2)) , 



a. 

a 



We focus on Jf^i^^ as the other cases are strictly analogous, indeed shghtly simpler. We 
have 

J- = 0,(1) 



^2Un 



and 



< 

I " 'I " \ IM / 

l=h 



E ^^Mt(^i+^2+/3+/4))i^lnNi'«'(f)-nNi'«'(#) 

hhlsk^'^'^'^^ l=h l=h 



< c 



N/8<h 



l=h l=h 
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Now notice that 

X1X2X3X4 - ymysy^ = X1X2X3 (0:4 - 1/4) + X1X2 {x3 - 1/3)1/4 + xi {x2 - 1/2)2/31/4 + {xi - 1/1)1/22/31/4 
whence 



max < —s-E 
N/s<h,i2,hM<N/2 y a% 



n \wi?a^{ 



< max 

N/fi<h,l2,hM<N/2 



+ max 

N/S<h,l2,l3M<N/2 



n 



i=h 



,2/4^' 



< 



i=h 



l~ |2 I |2 '3 i~i2 

FZil - Fur TT Fzi 



n 

i=h 



a 



N 



+ 



n 

i=h 



I~l2 I I2l 12 '2 |~|2 

F«3r - F/al Fur -Q Fzr 



'AT 



'N 



l=h 



'N 



+ 



+ max <; I Yla\%)]E 



'AT 



'AT 



'AT 



+ max 

N/%<h,h,h,k<N/2 



k 

n 



'N 



'N 



+ 



< 



c 

< —7- max 

Cat N/%<h,l2,h,k<N/2 



{E\wi\'' + E\wi\^Y X ^ 



max 



<J}q N/8<h,l2,l3,k<N/2 
0(1). 



□ 



Remark 7.4 It is evident from the proof that, in order to estabhsh (7.34), it is sufficient 
to impose the minimal bandwidth condition 

1 N 

— + — ^OasM^oo. 

M 

which is weaker than Assumption D. This is intuitively due to peculiar form that the 
aliasing effect assumes for Fourier transforms in the discrete case, see (7.32). 



8 Monte Carlo evidence 



In a way of confirmation of the CLT, simulations has been performed. It has been chose 
to put ci = /-I 1600 fields were simulated on the torus with these specifications and 
for each of them the skewness and kurtosis statistics were computed for = 2^^. The 
values obtained were then normalized dividing the by the theoretical standard deviations 
of these statistics, computed using formulas (3.12), (4.15), (4.16) and (4.17). Histograms 
are presented in Figures 2 and 3, showing a good accordance with Gaussianity. 



REFERENCES 



32 




-4-3-2-10 1 2 3 4 

Figure 2: Histogram of tiie skewness statistics over 1600 simulated fields. Here N = 2 




Figure 3: Histogram of the kurtosis statistics over the same 1600 simulated fields. 
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